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A correlated inhomogeneous mean-field approach is proposed in order to study a tight-binding 
model of the manganite heterostructures {LaMn03)2n/ {SrMn03)n with average hole doping x = 
1/3. Phase diagrams, spectral and optical properties of large heterostructures (up to 48 sites along 
the growth direction) with a single interface are discussed analyzing the effects of electron-lattice 
anti-adiabatic fluctuations and strain. The formation of a metallic ferromagnetic interface is quite 
robust with varying the strength of electron-lattice coupling and strain, though the size of the 
interface region is strongly dependent on these interactions. The density of states never vanishes 
at the chemical potential due to the formation of the interface, but it shows a rapid suppression 
with increasing the electron-lattice coupling. The in-plane and out-of-plane optical conductivities 
show sharp differences since the in-plane response has metallic features, while the out-of-plane one 
is characterized by a transfer of spectral weight to high frequency. The in-plane response mainly 
comes from the region between the two insulating blocks, so that it provides a clear signature of the 
formation of the metallic ferromagnetic interface. 

PACS numbers: 



I. INTRODUCTION 

Transition metal oxides are of great current inter- 
est because of the wide variety of the ordered phases 
that they exhibit and the strong sensitivity to external 
perturbationsii Among them, manganese oxides with for- 
mula Ri-x^xMnOs {R stands for a rare earth as La, A 
represents a divalent alkali element such as Sr or C a and 
X the hole doping), known as manganites, have been stud- 
ied intensively both for their very rich phase diagram and 
for the phenomenon of colossal magnetoresistanceiS This 
effect is often exhibited in the doping regime 0.2 < a: < 
0.5, where the ground state of the systems is ferromag- 
netic. The ferromagnetic phase is usually explained by 
invoking the double exchange mechanism in which hop- 
ping of an outer-shell electron from a Mn^'^ to a Mn'^'^ 
site is favored by a parallel alignment of the core spins <^ 
In addition to the double-exchange term that promotes 
hopping of the carriers, a strong interaction between elec- 
trons and lattice distortions plays a non-negligible role 
in these compounds giving rise to formation of polaron 
quasi-particles.— 

Very recently, high quality atomic-scale "digital" het- 
erostructures consisting of combination of transition 
metal oxide materials have been realized. Indeed, het- 
erostructures represent the first steps to use correlated 
oxide systems in realistic devices. Moreover, at the inter- 
face, the electronic properties can be drastically changed 
in comparison with those of the bulk. Recent examples 
include the formation of a thin metallic layer at the inter- 
face between band and Mott insulators as, for example, 
between SrTiO^ (STO) and LaTiOs oxides^ or between 
the band insulators^ LaAlOs and STO. 

Very interesting examples of heterostructure are given 



by the superlattices {LaMn03)m/ {SrMn03)n with 
n/{m + n) average hole doping.^ Here LaMnOa (LMO) 
(one electron per Mn Cg state) and SrMnO^ [SMO) 
(no electrons per Mn eg state) are the two end-member 
compounds of the alloy Lai^xSr^MnOs and are both 
antiferromagnetic insulating. In these systems, not only 
the chemical composition but also the thickness of the 
constituent blocks specified by m and n is important 
for influencing the properties of superlattices. Focus has 
been on the case m — 2n corresponding to the aver- 
age optimal hole doping x = 1/3^ The superlattices 
exhibit a metal-insulator transition as function of tem- 
perature for n < 2 and behave as insulators for n > 3. 
The superlattices undergo a rich variety of transitions 
among metal, Mott variable range hopping insulator, 
interaction-induced Efros-Shklovskii insulator, and pola- 
ronic insulator. ^'^ 

Interfaces play a fundamental role in tuning the metal- 
insulator transitions since they control the effective dop- 
ing of the different layers. Even when the system is 
globally insulating (n > 3), some nonlinear optical mea- 
surements suggest that, for a single interface, ferromag- 
netism due to double-exchange mechanism can be in- 
duced between the two antiferromagnetic blocks .ii More- 
over, it has been found that the interface density of 
states exhibits a pronounced peak at the Fermi level 
whose intensity correlates with the conductivity and 
magnetization. -'^^ These measurements point toward the 
possibility of a two-dimensional half-metallic gas for the 
double-layer— whose properties have been studied by us- 
ing ab-initio density functional approaches. However, 
up to now, this interesting two-dimensional gas has not 
been experimentally assessed in a direct way by using lat- 
eral contacts on the region between the LM O and SMO 
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blocks. 

In analogy with thin films, strain is another important 
quantity in order to tune the properties of manganite 
heterostructures. For example, far from interfaces, inside 
LMO, electron localization and local strain favor antifer- 
romagnetism and Cg {3z^ — r^) orbital occupation.— The 
magnetic phase in LMO is compatible with the C typei^ 
Moreover, by changing the substrate, the ferromagnetism 
in the superlattice can be stabilized.— 

From the theoretical point of view, in addition to ab 
initio calculations, tight-binding models have been used 
to study manganite super lattices. Effects of magnetic 
and electron-lattice interactions on the electronic proper- 
ties have been investigated going beyond adiabatic mean- 
field approximations . ^ ''i However, the double layer with 
large blocks of LMO and SMO has not been much stud- 
ied. Moreover, the effects of strain have been analyzed 
only within mean-field approachesi^ 

In this paper we have studied phase diagrams, spec- 
tral and optical properties for a very large bilayer 
{LM0)2nl {SMO)n (up to sizc of 48 planes relevant for 
a comparison with fabricated heterostructures) starting 
from a tight binding model. We have developed a cor- 
related inhomogeneous mean-field approach taking into 
account the effects of electron-lattice anti-adiabatic fluc- 
tuations. Strain is simulated by modulating hopping 
and spin-spin interaction terms. We have found that a 
metallic ferromagnetic interface forms for a large range of 
the electron-lattice couplings and strain strengths. For 
this regime of parameters, the interactions are able to 
change the size of the interface region. We find the 
magnetic solutions that are stable at low temperature 
in the entire superlattice. The general structure of our 
solutions is characterized by three phases running along 
growth z-direction: antiferromagnetic phase with local- 
ized/delocalized (depending on the model parameters) 
charge carriers inside LMO block, ferromagnetic state at 
the interface with itinerant carriers, localized polaronic 
G-type antiferromagnetic phase inside SMO block. The 
type of antiferromagnetic order inside LMO depends on 
the strain induced by the substrate. 

We have discussed the spectral and optical properties 
corresponding to different parameter regimes. Due to 
the formation of the metallic interface, the density of 
states is finite at the chemical potential. With increas- 
ing the electron-phonon interaction, it gets reduced at 
the chemical potential, but it never vanishes even in the 
intermediate to strong electron-phonon coupling regime. 
Finally, we have studied both the in-plane and out-of- 
plane optical conductivities pointing out that they are 
characterized by marked differences: the former shows a 
metallic behavior, the latter a transfer of spectral weight 
at high frequency due to the effects of the electrostatic 
potential well trapping electrons in LMO block. The 
in-plane response at low frequency is mainly due to the 
region between the two insulating blocks, so that it can 
be used as a tool to assess the formation of the metallic 
ferromagnetic interface. 



The paper is organized as follows: in sec. II the model 
and variational approach are introduced, in III the results 
regarding the phase diagrams are discussed, in sec. IV 
the spectral properties and in sec. V the optical conduc- 
tivities are analyzed, in the final section the conclusions. 

II. THE VARIATIONAL APPROACH 
A. Model Hamiltonian 

For manganite superlattices, the hamiltonian of the 
bulk Hf) has to be supplemented by Coulomb terms rep- 
resenting the potential arising from the pattern of the La 
and Sr ions,^^ thus 

H = Hf) + Hcoui- (1) 

In order to set up an appropriate model for the double 
layer, it is important to take into account the effects of 
the strain. The epitaxial strain produces the tetrago- 
nal distortion of the MuOq octahedron, splitting the eg 
states into — and — statesi^ If the strain is 
tensile, — is lower in energy, while, if the strain is 
compressive, "iz^ — is favored. In the case of n = 8 and 
three interfaces jii^ the superlattices grown on STO are 
found to be coherently strained: all of them are forced 
to the in-plane lattice parameter of substrate and to an 
average out-of-plane parameter c ~ 3.87A.^ - As a conse- 
quence, one can infer that LMO blocks are subjected to 
compressive strain (—2.2%) and SMO blocks to tensile 
strain (-1-2.6%). For the case of LMO block, the result- 
ing higher occupancy of 3^^ — enhances the out-of- 
plane ferromagnetic interaction owing to the larger elec- 
tron hopping out-of-plane. For the case of SMO block, 
the reverse occurs. A suitable model for the bilayer has to 
describe the dynamics of the eg electrons which in LMO 
block and SMO block preferentially occupy the more 
anisotropic 3z^ — orbitals and more isotropic — 
orbitals, respectively. For this reason, in this paper we 
adopt an effective single orbital approximation for the 
bulk manganite. 

The model for the bulk takes into account the double- 
exchange mechanism, the coupling to the lattice distor- 
tions and the super-exchange interaction between neigh- 
boring localized t2g electrons on Mn ions. The coupling 
to longitudinal optical phonons arises from the Jahn- 
Teller effect that splits the eg double degeneracy. Then, 
the Hamiltonian Hq reads: 

^0 = - Vi — ^^n — 

Ri,S \ / 
Ri Ri 

+\Y.^\s\^r. ■ ^r.+s- t'Y.''n.''R.- (2) 

Ri,S Ri 
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Here t^g^ is the transfer integral of electrons occupying 

eg orbitals between nearest neighbor (nn) sites, S^"^'^^ 
is the total spin of the subsystem consisting of two local- 
ized spins on nn sites and the conduction electron, S*^ 

is the spin of the t2g core states (5 = 3/2), ct (cj^^ 

creates (destroys) an electron with spin parallel to the 
ionic spin at the i-th site in the eg orbital. The coordi- 
nation vector 6 connects nn sites. The first term of the 
Hamiltonian describes the double-exchange mechanism 
in the limit where the intra-atomic exchange integral J 
is far larger than the transfer integral t^g^. Furthermore, 
in cq.Q, uq denotes the frequency of the local optical 
phonon mode, aj^ ('^^ ) creation (annihilation) 

phonon operator at the site i, the dimensionless param- 
eter g indicates the strength of the electron-phonon in- 
teraction. Finally, in Eq.(l2|), e|^-j represents the antifer- 
romagnetic super-exchange coupling between two nn t2g 
spins and fi is the chemical potential. The hopping of 
electrons is supposed to take place between the equivalent 
nn sites of a simple cubic lattice (with finite size along the 
z axis corresponding to the growth direction of the het- 
erostructure) separated by the distance |n — n'| = a. The 
units are such that the Planck constant the Boltz- 

mann constant kB=^ and the lattice parameter a=\. 

Regarding the terms due to the interfaces, one con- 
siders that La?^ and Sr^^ ions act as -1-1 charges of 
magnitude e and neutral points, respectively. In the het- 
erostructure, the distribution of those cations induces an 
interaction term for eg electrons of Mn giving rise to the 
Hamiltonian 



a = 4A, and typical value of dielectric constant e ^ 10 
to be around 0.2. 

Strain plays an important role also by renormalizing 
the heterostructure parameters. Strain effects can be 
simulated by introducing an anisotropy into the model 
between the in-plane hopping amplitude — t (with 
5\\ indicating nearest neighbors in the x — y planes) and 
out-of-plane hopping amplitude t|5^| = tz (with 5z indi- 
cating nearest neighbors along z axis)i^ Moreover, the 
strain induced by the substrate can directly affect the 
patterns of core spins. Therefore, in our model, we 
have also considered the anisotropy between the in-plane 
super-exchange energy ei^n | = e and the out-of-plane one 
— £2- We have found that the stability of magnetic 
phases in LMO blocks is influenced by the the presence 
of compressive strain, while in SMO the sensitivity to 
strain is poor. Therefore, in all the paper, we take as 
reference the model parameters of the SMO layers and 
we will consider anisotropy only in the LM O blocks with 
values of the ratio tz/t larger than unity and of the ratio 
e^/e smaller than unity. 

Finally, in order to investigate the effects of the 
electron-lattice coupling, we will use the dimensionless 
quantity A defined as 
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with n^^ ~ '^R electron occupation number at Mn 

site i, Ri and .Rf" are the positions of Mn and La'^+ 
in ith unit cell, respectively, and Cd is the dielectric con- 
stant of the material. In our calculation the long-range 
Coulomb potential has been modulated by a factor rj in- 
ducing a fictitious finite screening-length (see Appendix 
A). This factor was added only for computational rea- 
sons since it allows to calculate the summations of the 
Coulomb terms over the lattice indices. We have mod- 
eled the heterostructures as slabs whose in-plane size is 
infinite. 

In order to describe the magnitude of the Coulomb 
interaction, we define the dimensionless parameter a = 
e^/(ae£;i|j|) which controls the charge-density distribu- 
tion. The order of magnitude of a can be estimated from 
the hopping parameter t^g^ ~ 0.65eV, lattice constant 
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(3) 



A 



(4) 



In all the paper we will assume oJo/t — 0.5. 



B. Test Hamiltonian 

In this work, we will consider solutions of the hamilto- 
nian that break the translational invariance in the out- 
of-plane z-direction. The thickness of the slab is a pa- 
rameter of the system that will be indicated by Nz- We 
will build up a variational procedure including these fea- 
tures of the heterostructures. A simplified variational 
approach similar to that developed in this work has al- 
ready been proposed by some of the authors for mangan- 
ite bulks^^ and filmsi^i^ 

In order to treat variationally the electron-phonon in- 
teraction, the Hamiltonian ([T]) has been subjected to an 
inhomogeneous Lang-Firsov canonical transformation.^— 
It is defined by parameters depending on plane indices 
along z-direction: 



U = exp 



E(/- 



. . (5) 

where i|| indicates the in-plane lattice sites (ix,iy), 
while iz the sites along the direction z. The quantity 
fi^ represents the strength of the coupling between an 
electron and the phonon displacement on the same site 
belonging to i^-plane, hence it measures the degree of the 
polaronic effect. On the other hand, the parameter A^^ 
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denotes a displacement field describing static distortions 
that are not influenced by instantaneous position of the 
electrons. 

In order to obtain an upper limit for free energy, the 
Bogoliubov inequality has been adopted: 



F < Ftest + {H — Htest)t, 



(6) 



where Ff^st a-nd Htf,gt are the free energy and the Hamil- 
tonian corresponding to the test model that is assumed 
with an ansatz. H stands for the transformed Hamilto- 
nian H — UHW . The symbol {)t indicates a thermody- 
namic average performed by using the test Hamiltonian. 
The only part of Htest which contributes to {H — Htest)t 
is given by the spin freedom degrees and depends on the 
magnetic order of the t2g core spins. For the spins, this 
procedure is equivalent to the standard mean-field ap- 
proach. 

The model test hamiltonian, Htest, is such that that 
electron, phonon and spin degrees of freedom are not 
interacting with each other: 



(7) 



The phonon part of Htest simply reads 

Hfest = ^0 I I 'i'^ ' 

and the spin term is given by 

HZt = E E I I ' (9) 



where gs is the dimensionless electron-spin factor {gs — 
2), iJB is the Bohr magneton, and hf^^^^^ is the effective 
variational magnetic field. In this work, we consider the 
following magnetic orders modulated plane by plane: 

F, h^^^^^Kl; 

G, /if,,,,, = (-ir+^^+^i/iti. (10) 

For all these magnetic orders, the thermal averages of 
double-exchange operator, corresponding to neighboring 
sites in the same plane iz 7i,;i||,i|j-i-5|| and in different 
planes r]i^,iz+s^;i^^, preserve only the dependence on the 
z plane index: 

Tiilill ,i| |-|-i5|| ^ \ r)C I 1 /* ~ Tii 



Tl — / 



25+1 
1/2 



2S+1 



)t = V^.,^.+S,■ (11) 



In order to get the mean-field electronic Hamiltonian, 
we make the Hartree approximation for the Coulomb in- 
teraction. The electronic contribution H^lg^ to the test 
Hamiltonian becomes 



m. 



»|| i.=l <5| 



EE i.e., , 



ill i^ = l 



(12) 

In Eq. p^ . the quantity (peff{iz) indicates the effective 
potential seen by the electrons. It consists of the Hartee 
self-consistent potential 0(iz) (see Appendix A) and a 
potential due to the electron-phonon coupling: 



4>eff{iz) = (l>{iz) +5^woC,,, (13) 



with 



a. 



/2 -2/,, +2A,,(/,, -1). (14) 



The factors 



and 



represent the phonon 



thermal average of Lang-Firsov operators: 

=(^.|,^.^/,„..+^.)t, (15) 
where the operator Xj^ reads 

Qfi^ (fi)3 — a^- ) 

R^ ~ 

Finally, the quantity Ti and T2 derive from the Hartree 
approximation (see Appendix A), and Ny denote the 
size of the system along the two in-plane directions, re- 
spectively. In order to calculate the variational free en- 
ergy, we need to know eigenvalues and eigenvectors of 
H^l^^ which depend on the magnetic order of core spins 
through the double exchange terms. 



Magnetic order and diagonalization of the 
electronic mean-field Hamiltonian 



In order to develop the calculation, we need to fix the 
magnetic order of core spins. The patters of magnetic or- 
ders is determined by the minimization of the total free 
energy. By exploiting the translational invariance along 
the directions perpendicular to the growth axis of the 
heterostructure, the diagonalization for H^l.g^ reduces to 



an effective unidimensional problem for each pair of con- 
tinuous wave vectors {kx,ky) = k^^. For some magnetic 
patterns, the electronic problem is characterized at the 
interface by a staggered structure. Therefore, we study 
the electron system considering a reduced first Brillouin 
zone of in-plane wave vectors. To this aim, we represent 
ff^g'^j with the 2N^ states 



\kxjky,iz), \kx -\- , ky -\- TT , iz) , 



(16) 



with the wave vectors such that — 7r/2 < k^ < 7r/2, 
— 7r/2 < ky < 7r/2, and iz going from 1 to Nz- The 
eigenstates of electronic test Hamiltonian are indicated 
by E{kx,ky,n), with the eigenvalue index n going from 
1 to 2Nz ■ The eigenvector related to n is specified in the 
following way: 6i^(fc||,n) for the first Nz components, 

Pi^{k\\,n) for the remaining Nz components. 

The variational procedure is self-consistently per- 
formed by imposing that the total density of the system 
p is given by N^a/Nz, with N^a the number of layers of 
LMO block, and the local plane density xi^z) is equal to 
(n^,). Therefore, one has to solve the following iV^ -|- 1 
equations: 



1 



NxNyN, 



(17) 



and 



NxNy - 

k 



£;(A?||,n) 



\b,^{ku,n)\^ + \pi^{ku,n)\^ + 



[bt I , ")k, (fc| I , n) -f p*Jfc| I , n)&i^ (A;| I , n)] 



(18) 

where np{z) is the Fermi distribution function. These 
equations allow to obtain the chemical potential fi and 
the the local charge density xi'i'z)- As result of the varia- 
tional analysis, one is able to get the charge density pro- 
file corresponding to magnetic solutions which minimize 
the free energy. 



III. STATIC PROPERTIES AND PHASE 
DIAGRAMS 



We have found the magnetic solutions and the corre- 
sponding density profiles that are stable for different sizes 
of the LMO and SMO blocks. The inhomogeneous vari- 
ational approach allows to determine the values of the 
electron-phonon parameters /i. and A^. , and the mag- 
netic order of the t2g spins through the effective magnetic 
fields hi^ . We will study the systems in the interme- 
diate to strong electron-phonon regime characteristic of 
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FIG. 1: Comparison among density profiles corresponding to 
different sizes at A = 0.5 and e = O.Olt. The index indicates 
the interface Afn-plane between the last La-plane in LMO 
block and the first Sr-plane in SMO block. 



manganite materials focusing on two values of coupling: 
A = 0.5 and A ~ 0.8. The maximum value of in-plane 
antiferromagnetic super-exchange is e = O.Olt. The value 
of the Coulomb term a is fixed to a = 0.2. We will ana- 
lyze the heterostrucures in the low-temperature regime: 
T = O.OSt. 

The general structure of our solutions is characterized 
by three phases running along z-direction. Actually, ac- 
cording to the parameters of the model, we find G or 
C antiferromagnetic phases corresponding to localized 
or delocalized charge carriers inside LMO block, respec- 
tively. The localization is ascribed to the electron-phonon 
coupling which gives rise to the formation of small po- 
larons. For the values of A considered in this paper, a fer- 
romagnetic phase always stabilizes around the interface. 
The size of the ferromagnetic region at the interface is 
determined by the minimization of the free energy and 
depends on the values of the system parameters. Only 
for larger values of A and e, the possibility of interface 
ferromagnetism is forbidden. Inside the SMO block, a 
localized polaronic G-type antiferromagnet phase is al- 
ways stable. 

At first, we have analyzed the scaling of the static 
properties as function of the size of the system along 
the z growth direction. Therefore, a comparison of the 
density profiles has been done with {LMO)%/ {SMO)i^ 
{LMO)ie/{SMO)8 and {LMO)32/{SMO)i6 systems. In 
Fig. [1] we show the density profiles in a situation where 
strain-induced anisotropy has not been introduced. It is 
worth noticing that we indicate the interface Mn-plane 
between the last La-plane in LMO block and the first Sr- 
plane in SMO block with the index 0. For a sufficiently 
large numbers of planes, the charge profile along z shows 
a well-defined shape. Indeed, the local density is nearly 



unity in LMO block, nearly zero in SMO block, and it 
decreases from 1 to in the interface region. The decrease 
of charge density for the first planes of LMO is due to the 
effect of open boundary conditions along the z direction. 
In the intermediate electron-phonon coupling regime that 
we consider in Fig. ([T]), the region with charge dropping 
involves 4 — 5 planes between the two blocks. We notice 
that the local charge density for {LMO)iqI {SMO)% and 
(LAfO)32/(S'AfO)i6 systems are very similar around the 
interface. Furthermore, the numerical results show close 
values of variational free energy corresponding to above 
mentioned systems. Given the similarity of the proper- 
ties of these two systems, in the following, we will develop 
the analysis on the role of interface studying the system 
[LMO)w/{SMO)8. 

For the same set of electron-phonon and magnetic cou- 
plings, the variational parameters and the Hartree self- 
consistent potential along z-axis are shown in Fig. 2. The 
effective magnetic fields are plotted for the most stable 
magnetic solution: antiferro G orders well inside LMO 
(planes 1 — 15) and SMO (planes 19 — 24), and ferromag- 
netic planes at the interface (planes 16 — 18). The peak 
in the plot of the magnetic fields signals that ferromag- 
netism is quite robust at the interface. The variational 
electron-phonon parameters fi^ are small on the LMO 
side and at the interface, but close to unity in SMO 
block. This means that, for these values of the couplings, 
carriers are delocalized in LMO up to the interface re- 
gion, but small polarons are present in the SM O block. 
The quantities A^^ , entering the variational treatment 
of the electron-phonon coupling, are determined by fi^ 
and the local density < rii^ > through the equation: 
Ai^ =< rii^ > (1 — /i,). The Hartree self-consistent 
potential $ indicates that charges are trapped into a po- 
tential well corresponding to the LMO block. Moreover, 
it is important to stress the energy scales involved in the 
well: the barrier between LMO and SMO block is of 
the order of the electron band-width. Furthermore, at 
the interface, the energy difference between neighboring 
planes is of the order of the hopping energy t. 

As mentioned above, for these systems, strain plays an 
important role. In order to study quantitatively its ef- 
fect, we have investigated the phase diagram under the 
variation of the hopping anisotropy tz/t for two differ- 
ent values of {ez ~ e — O.Olt, ez — 0). Indeed, we 
simulate the compressive strain in the LMO block in- 
creasing the ratio tz/t and decreasing Cz/e. On the other 
hand, the tensile strain in the SMO block favour the 
more isotropic — orbital and does not yield sizable 
effects. Therefore, for the SMO block, in the following, 
we choose tz = t and = e. For what concerns the 
electron-phonon interaction, we assume an intermediate 
coupling, A = 0.8. As shown in the upper panel of Fig. 
3, with increasing the ratio tz/t up to 1.7 for = e, 
the magnetic order in LMO does not change since it re- 
mains G antiferromagnetic. However, the character of 
charge carriers is varied. Actually, for A = 0.8, in the 
absence of anisotropy, small polarons are present in the 
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FIG. 3: Phase diagram in the hopping anisotropy-energy 
plane for LMOieSMOg system, corresponding to A = 0.8 
for — O.Oli (upper panel) and — (lower panel). 



LMO block. Moreover, at tz/t ~ 1.5, in LMO, a change 
from small localized polarons to large delocalized polaron 
occurs. For all values of the ratio tz/t, the interface re- 
gion is characterized by ferromagnetic order with large 
polaron carriers and SMO by G antiferromangnetic or- 
der with small polaron carriers. 

It has been shown that it is also important to consider 
the anisotropy in super-exchange [tz 7^ e) parameters 
as consequence of strain.— In order to simulate the ef- 
fect of compressive strain in LMO, a reduction of will 
be considered. We discuss the limiting case: €z = 0. 
For this regime of parameters, the effect on the magnetic 



phases is the strongest. As shown in the lower panel of 
Fig. 3, for 1.28 < tjt < 1.5, in LMO block, a C-type 
antiferromagnetic phase is the most favorable. The tran- 
sition from small to large polaron again takes place at 
tz/t ~ 1.5. Therefore, we have shown that there is a 
range of parameters where LMO block has C-type anti- 
ferromagnetic order with small localized polarons. Due 
to the effect of strain, the magnetic solution in LMO 
turns out to be compatible with experimental results in 
superlatticesji^ The interface is still ferromagnetic with 
metallic large polaron features. In the figure A/B/C 
refers to magnetic orders and character of charge carriers 
inside LMO (A), at interface (B), inside SMO (C). 

In order to analyze the effects of the electron-phonon 
interaction, a comparison between two different electron- 
phonon couplings is reported in Fig. 4. We have inves- 
tigated the solutions which minimize the variational free 
energy at fixed value of the anisotropy factors tz/t = 1.3 
and Cz = at X = 0.5 and A — 0.8. The magnetic solu- 
tion in LMO block is C antiferromagnetic until the 15th 
plane. For both values of A, polarons are small. In SMO 
block, starting from the 19th plane, the solution is G- 
type antiferromagnetic together with localized polarons. 
Three planes around the interface are ferromagnetically 
ordered. For A = 0.5, all the three planes at the inter- 
face are characterized by delocalized polarons, while, for 
A = 0.8, only the plane linking the ends of LMO and 
SMO blocks is with delocalized charge carriers. 

As shown in Fig. 4, the quantity A has important 
consequences on the physical properties such as the lo- 
cal particle density. Actually, for A = 0.8 the transition 
from occupied to empty planes is sharper at the interface. 
Only one plane at the interface shows an intermediate 
density close to 0.5. For A = 0.5 the charge profile is 
smoother and the three ferromagnetic planes with large 
polarons have densities different from zero and one. 

For the analysis of the spectral and optical quantities, 
we will consider the parameters used for the discussion 
of the results in this last figure. 



IV. SPECTRAL PROPERTIES 

In the following section we will calculate the spectral 
properties of the heterostructure for the same parameters 
used in Fig. 4. 

Performing the canonical transformation ([SJ and ex- 
ploiting the cyclic properties of the trace, the electron 
Matsubara Green's function becomes 



G{Ri, Rj 



-(T,C5.(r)X. (r)4 (0)xt (0)). (19) 



By using the test Hamiltonian (O, the correlation func- 
tion can be disentangled into electronic and phononic 
termsi^i^ Going to Matsubara frequencies and making 
the analytic continuation za;„ tu + id, one obtains the 
retarded Green's function and the diagonal spectral func- 
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FIG. 4: Comparison between local particle density corre- 
sponding to A = 0.5 and A = 0.8. 



tion A^^^^ [ui) corresponding to Ri — Rj 
Ai^/^uj) = 

OO 

e^T Ii{S'^)e^[l~nF{uj-liOo)]gly''{io-liOo) 



1— — 00 



I— — CO 

(20) 

where = g^fl{2No + 1), S'^^ = 2g^fl [No{No + 1)]^ , 
Ii{z) modified Bessel functions, and (7*'"'*"(aj) is 



27T 



2N^ 



|6i^(fc||,n)p + |p,Jfc||,n)p + 

(21) 

The density of states D{uj) is defined as 
1 1 



D{cj) 



N^NyNz 27r 



(22) 



In Fig. 5 we report the density of state of the system 
{LMO)i(,/ {SMO)s- It has been calculated measuring 
the energy to the chemical potential /i. This comparison 
has been made at fixed low temperature [KbT — 0.05t), 
therefore we can consider the chemical potential very 
close to the Fermi energy of the system. At A = 0.5, 
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even a weak disorder induces localization. On the other 
hand, a weak disorder is not able to prevent the forma- 
tion of the ferromagnetic metallic interface favored by the 
double-exchange mechanism and the charge transfer be- 
tween the bulk-like blocks: the states at the Fermi level 
due to the interface formation have enough density^'^ so 
that they cannot be easily localized by weak disorder. In 
this section, we have shown that this can be the case in 
the intermediate electron-phonon coupling regime appro- 
priate for LMO/SMO heterostructures. 

In the next section we will analyze the effects of 
electron-phonon coupling and strain on the optical con- 
ductivity in the same regime of the parameters considered 
in this section. 



FIG. 5: Comparison between density of states (in units of 
1/t) as a function of the energy (in units of t) corresponding 
to A = 0.5 and A = 0.8. 



the spectral function exhibits a residual spectral weight 
at jU. The main contribution to the density of states at 
the chemical potential fi comes from the three ferromag- 
netic large polaron planes at the interface. Indeed, the 
contributions due to the (LMO) and (SMO) blocks is 
negligible. 

For stronger electron-phonon coupling at A = 0.8, we 
observe an important depression of the spectral function 
at ^. Hence the formation of a clear pseudogap takes 
place. This result is still compatible with the solution of 
our variational calculation since, for this value of A = 0.8, 
there is only one plane with delocalized charge carriers 
which corresponds to the plane indicated as the interface 
{iz — 17), while the two further ferromagnetic planes 
around the interface are characterized by small polarons. 
The depression of the density of the states at the Fermi 
energy is due also to the polaronic localization well inside 
the LMO and SM O block. In any case we find that, even 
for A = 0.8, the density of states never vanishes at the 
interface in agreement with experimental results 

In this section we have found strong indications that a 
metallic ferromagnetic interface can form at the interface 
between LMO and SMO blocks. This situation should 
be relevant for superlattices with n > 3, where resistiv- 
ity measurements made with contacts on top of LMO 
show a globally insulating behavior. In our analysis we 
have completely neglected any effect due to disorder even 
if, both from experiments^!^ and theories^i, it has been 
suggested that localization induced by disorder could be 
the cause of the metal-insulator transition observed for 
n > 3. We point out that the sizable source of disorder 
due to the random doping with S'r^+ is strongly reduced 
since, in superlattices, ia'^+ and S'r^+ ions are spatially 
separated by interfaces. Therefore, the amount of dis- 
order present in the heterostructure is strongly reduced 
in comparison with the alloy. However, considering the 
behavior of the LMO (SMO) block as that of a bulk 
with a small amount of holes (particles), one expects that 



V. OPTICAL PROPERTIES 

To determine the linear response to an external field 
of frequency w, we derive the conductivity tensor CTq^^ 
by means of the Kubo formula. In order to calculate the 
absorption, we need only the real part of the conductivity 



1 IILLL^ 



(23) 



where H^*^^ is the retarded current-current correlation 
function. Following a well defined scheme^^i^i and ne- 
glecting vertex corrections, one can get a compact ex- 
pression for the real part of the conductivity aa,a- It 
is possible to get the conductivity both along the plane 
perpendicular to growth and parallel to it, azz- 

In order to calculate the current-current correlation func- 
tion, one can use the spectral function Ar ■ ■ derived 
in the previous section exploiting the translational in- 
variance along in-plane direction. It is possible to show 
that the components of the real part of the conductivity 
become 



Re[axx]iuj) = 



~( ^^[nF(wi - W) - 71^(^1 )] 

w J_oo 47r 

j^{uJi -a;)ylfc^,fe„;,^jjwi), (24) 



X 



and 



E 1^ E E 

dioi 



^ kx ,ky -Iz ,Jz Oiz ,(); 



^ -l-oo 47r 



(25) 



In Fig. 6, we report the in-plane conductivity as func- 
tion of the frequency at A = 0.5 and A = 0.8. We have 
checked that the in-plane response mainly comes from 
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FIG. 6: The conductivity (in units of /{rat), with m — 
\/{2t)) into the plane perpendicular to growth direction of 
the {LMO)i6/{SMO)s bilayer as a function of the energy (in 
units of t) for different values of A. 

the interface planes. Both conductivities are character- 
ized by a Drude-like response at low frequency. There- 
fore, the in-plane conductivity provides a clear signature 
of the formation of the metallic ferromagnetic interface. 
However, due to the effect of the interactions, we have 
found that the low frequency in-plane response is at least 
one order of magnitude smaller than that of free electrons 
in the heterostructures. Moreover, additional structures 
are present in the absorption with increasing energy. For 
A = 0.5, a new band with a peak energy of the order of 
hopping t = 2wo is clear in the spectra. This structure 
can be surely ascribed to the presence of large polarons 
at the three interface planes.™ Actually, this band comes 
from the incoherent multiphonon absorption of large po- 
larons at the interface. This is also confirmed by the fact 
that this band is quite broad, therefore it can be inter- 
preted in terms of multiple excitations. For A = 0.8, the 
band is even larger and shifted at higher energies. In 
this case, at the interface, large and small polarons are 
present with a ferromagnetic spin order. Therefore, there 
is a mixing of excitations whose net effect is the transfer 
of spectral weight at higher frequencies. 

The out-of-plane optical conductivities show signif- 
icant differences in comparison with the in-plane re- 
sponses. In Fig. 7, we report out-of-plane conductivity as 
function of the frequency at A = 0.5 and A = 0.8. First, 
we observe the absence of the Drude term. Moreover, the 
band at energy about 2ujq is narrower than that in the 
in-plane response. Therefore, the origin of this band has 
to be different. Actually, the out-of-plane optical con- 
ductivities are sensitive to the interface region. A charge 
carrier at the interface has to overcome an energy barrier 
in order to hop to the neighbour empty site. As shown 
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FIG. 7: The conductivity (in units of e^/(mt), with m — 
l/(2t)) along the growth direction of the {LMO)i6{SMO)s 
bilayer as a function of the energy (in units of t) for A — 0.5 
and A = 0.8. 



in Fig. 2, the typical energy for close planes at the inter- 
face is of the order of the hopping t. Therefore, when one 
electrons hops along z, it has to pay at least an energy 
of the order of t. In the out-of-plane spectra, the peaks 
at low energy can be ascribed to this process. Of course, 
by paying a larger energy, the electron can hop to next 
nearest neighbors. This explains the width of this band 
due to inter-plane hopping. 

Additional structures are present at higher energies in 
the out-of-plane conductivities. For A = 0.5 the band 
at high energy is broad with small spectral weight. For 
A — 0.8, there is an actual transfer of spectral weight at 
higher energies. A clear band is peaked around lOt. This 
energy scale can be intepreted as given by 2g^uJo = 9.6t 
for A — 0.8. Therefore, in the out-of-plane response, the 
contribution at high energy can be interpreted as due to 
small polarons. -^i^I 

Unfortunately, experimental data about optical prop- 
erties of the LMO / SMO bilayers are still not available. 
Therefore, comparison with experiments is not possible. 
Predictions about the different behaviors among Uxx and 
Gzz can be easily checked if one uses in-plane and out- 
of-plane polarization of the electrical fields used in the 
experimental probes. More important, the formation of 
two-dimensional gas at the interface expects to be con- 
firmed by experiments made by using lateral contacts di- 
rectly on the region between the LMO and SMO blocks. 
The d.c. conductivity of the sheet could directly measure 
the density of carriers of the interface metal and con- 
firm the Drude-like low frequency behavior of in-plane re- 
sponse. Finally, one expects that a weak disorder present 
in the system and not included in our analysis can in- 
crease the scattering rate of the carriers reducing the 
value of the in-plane conductivity for a; — >■ 0. 
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VI. CONCLUSIONS 



In this paper we have discussed phase diagrams, 
spectral and optical properties for a very large bilayer 
{LM0)2n / {SMO)n (up to 48 sites along the growth 
direction). A correlated inhomogeneous mean- field ap- 
proach has been developed in order to analyze the effects 
of electron-lattice anti-adiabatic fluctuations and strain. 
We have shown that a metallic ferromagnetic interface is 
a quite robust feature of these systems for a large range of 
the electron-lattice couplings and strain strengths. Fur- 
thermore, we have found that the size of the interface 
region depends on the strength of electron-phonon inter- 
actions. At low temperature, the general structure of 
our solutions is characterized by three phases running 
along growth 2:-direction: antiferromagnetic phase with 
localized/Delocalized charge carriers inside LMO block, 
ferromagnetic state with itinerant carriers at the inter- 
face, localized polaronic G-type antiferromagnetic phase 
inside SMO block. The type of antiferromagnetic or- 
der inside LMO depends on the strain induced by the 
substrate. 

Spectral and optical properties have been discussed for 
different parameter regimes. Due to the formation of 
the metallic interface, even in the intermediate to strong 
electron-phonon coupling regime, the density of states 
never vanishes at the chemical potential. Finally, in- 
plane and out-of-plane optical conductivities are sharply 
different: the former shows a metallic behavior, the lat- 
ter a transfer of spectral weight at high frequency due 
to the effects of the electrostatic potential well trapping 
electrons in LMO block. The in-plane response provides 
a signature of the formation of the metallic ferromagnetic 
interface. 

In this paper we have focused on static and dynamic 
properties at very low temperature. The approach used 
in the paper is valid at any temperature. Therefore, it 
could be very interesting to analyze not only single in- 
terfaces, but also superlattices with different unit cells at 
finite temperature. Work in this direction is in progress. 



N, 



iz i\\ is — 1 

The self-consistent Hartee potential is given by 

.=2 



■[E x{j.)S{i.-3.) + 



E x{jz)S{i,-j.) + Smxiiz) - S^iQ], (A2) 



jz <iz 

where the quantity Ti is 

2 Nz Nz 



Ti = ^[E E XizX3zS{iz- jz) 



Nz 



Nz 



E X^zXjzS{^z - Jz) + S^{Q)J2xl], (A3) 



and T2 



2 NLa Nl 



^2 = ^[E E s{h-Jz) + 



NLa 



J2 S{h-J,) + NLaSl] 



(A4) 



with S{nz), S'i(O) end S2{nz) obtained by adding the 
Coulomb terms on in-plane lattice index. The summa- 
tions have been made modulating the Coulomb inter- 

.,, . „ , „2 p2 ->7s|Pi-»',| 

action with a screening factor: .^^ ^ , — — , 

° \ri-rj\ |n-rj| ' 

where — is a fictitious finite screening length in units of 
the lattice parameter a. Therefore, S{nz) is 



Appendix A 



exp i-risJml + ml + nl) 
S{nz)= 1 2 2 2 ' ^^^^ 



In this Appendix we give some details about the ef- 
fective electronic Hamiltonian derived within our ap- 
proach. After the Hartree approximation for the long- 
range Coulomnb interactions, the mean-field electronic 
Hamiltonian reads: 



51(0) is given by 



exp (-vsJml+ml) 
5r(0)= J2 , . ' (A6) 
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with {rrixTmy) ^ (0,0), and S2{iz — jz) is References 

exp i-risJhl + hl + hi) 

with Iz number of the planes of LMO block, hx = rux — 
0.5, hy — iriy — 0.5, and hz = Uz — iz — 0.5. 
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